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In our recent works we derived a chiral 0(q 4 ) two-pion exchange nucleon-nucleon potential (TPEP) 
formulated in a relativistic baryon (RB) framework, expressed in terms of the so called low energy 
constants (LECs) and functions representing covariant loop integrations. We showed that the ex- 
pansion of these functions in powers of the inverse of the nucleon mass reproduces most of the terms 
of the TPEP derived from the heavy baryon (HB) formalism, but such an expansion is ill defined 
and does not converge at large distances. In the present work we perform a study of the phase shifts 
in nucleon-nucleon (NN) scattering for peripheric waves (L > 3), which are sensitive to the tail of 
the potential. We assess quantitatively the differences between the RB and HB results, as well as 
f*^ ■ variations due to different values of the LECs. By demanding consistency between the LECs used 

in ttN and NN scattering we show how NN peripheral phase shifts could constrain these values. 
C^l , We demonstrate that this idea, first proposed by the Nijmegen group, favors a smaller value for the 

LEC C3 than the existing ones, when considering the TPEP up to order q 4 in the chiral expansion. 
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INTRODUCTION 



One of the most intriguing mysteries in the history of modern physics is trying to describe and understand nuclear 
' forces in terms of a basic, elementary theory of strong interactions. Despite more than fifty years of intense research 
t^J- \ in this area only in the last decade such a connection could be established, between what is believed to be the theory 

■ of strong interactions, QCD, and the rich phenomenology of low-energy nuclear physics. Chiral symmetry plays an 
essencial role in this game. Weinberg first pointed out that it can be used to obtain not only the same predictions 
of current algebra , but also corrections to them. Later these ideas were systematically elaborated by Gasser and 
Leutwyler Q, evolving to what is known as chiral perturbation theory (ChPT). It relies on the fact that the QCD 
lagrangian exhibits chiral SU(Nf) x SU(Nf) symmetry, Nf being the number of light quarks (u, d, and, eventually, s 
quarks), in the limit where their masses, m q , go to zero. This underlying theory is then described effectively including 

I , the relevant degrees of freedom (pions and nucleons for the case we are interested in, Nf = 2) in the most generic way 
• constrained only by the allowed symmetries, and caracterized by an expansion parameter q representing either the 
^ | momenta involved or the symmetry breaking terms m q cx /i 2 , fi being the Goldstone boson (pion) mass. The effective 

■ theory is then evaluated by means of feynman diagrams, strictly obeying a given power counting scheme. These ideas 
' were successfuly applied to many processes in the mesonic sector, forming a compelling evidence that ChPT indeed 

• i— i . represents QCD in the low-energy regime 0,0,0- 

' The power counting is an essential ingredient in ChPT, allowing us to evaluate corrections in a consistent and 
^ \ controlled way. In the mesonic sector it follows naturally, once loop contributions are regularized using dimensional 

- - 1 regularization. On the other hand, with baryons an additional energy scale (the baryon mass) is introduced, which 
does not vanish in the chiral limit and makes the theory much more complicated, with chiral loops (using dimensional 
regularization) spoiling the power counting • For many years it was believed that a relativistic treatment of baryons 
in ChPT was not consistent with power counting. The first and widely used formalism which restores the latter 
consists in applying a sort of non-relativistic expansion at the level of the lagrangian, around the limit of infinitely 
heavy baryon (HBChPT), proposed by Jenkins and Manohar in the early 90s j8j. Almost ten years later, Becher 
and Leutwyler 0, Il0l | showed that it is possible to formulate baryonic ChPT preserving explicit Lorentz invariance 
and power counting (RBChPT). Based on previous works of Ellis and Tang [llj, the so called infrared regularization 
allows one to not only remove the power counting violating terms from loop integrals, but also recover the results 
from HBChPT when subjected to an expansion in powers of the inverse of the baryon mass. From the latter, Becher 
and Leutwyler showed that the heavy baryon approach fails to preserve the correct analytic structure of the theory 
in a specific low energy domain. In irN scattering, this region corresponds to the Mandelstam variable t close to 4/x 2 . 



This motivated further studies inside the covariant framework 



lorrcsponcis 

EH, HI IH 



* Email address: higa@jlab.org 



In recent works 0, ^| we brought the issue of HBChPT and RBChPT to the two-pion exchange nucleon-nucleon 
potential (TPEP), showing that its long distance properties are governed by the same low energy region where the 
HBChPT has convergence problems. In we also performed an expansion of our results in powers of 1/m (m 
being the nucleon mass) and compared with the HB expressions from Entem and Machleidt |17| . where we found a 
small number of differences. These works triggered the importance to assess quantitatively the differences between 
the HB- and RB-ChPT NN potentials in asymptotic phase shifts, and if they are comparable to the uncertainties 
in the parameters of the potential, the so called low energy constants (LECs). The discrepancies found between our 
expanded fl5l ] and the HB |l7j expressions also needs to be understood as they should, in principle, coincide. It is the 
main goal of this work to address both issues. 

The determination of the LECs is of much current interest, as they describe the dynamics of QCD at low energies. 
At the order we are working on, it requires the knowledge of the dimension two (c,-) and three (dj) LECs, present 

at the order q 2 and order q 3 pion-nuclcon chiral Lagrangians, respectively. Obviously, it seems more 

natural trying to obtain their values from ttN processes. In Ref. [lg it was performed through a fit to threshold and 
subthreshold parameters, and a physical interpretation for their values was given, similar to what has been done in the 
pure mesonic sector |l9j . For instance, assuming that the LEC ci is entirely saturated by a scalar meson exchange, a 
perfect agreement for their mass and coupling constant with the ones used by the Bonn one boson exchange model |20j 
emerges. The other dimension two LECs C2, C3, and C4 where also estimated to be dominated by the delta resonance, 
with relevant contributions from rho (04) and sigma (03) mesons, and marginal contributions from the Roper resonance 
7V*(1440). However, de scrip tion of irN phase shifts, albeit consistent, was not satisfactory, the same happening with 
the values from Moizis |2lL l22l|. It was improved in subsequent papers [23I |24| . but the central values for the LECs 
still vary considerably 1251 . Moreover, the most reasonable results were obtained using the experimental analyses from 
the Karlsruhe group |26l |27|. which show disagreements with most recent irN database. As an alternative to irN 
processes, recently the Nijmegen group proposed the extraction of the dimension two LECs c\, C3, and C4 using NN 
partial wave analyses |28l |29| , where they obtained values with (statistical) uncertainties much smaller than obtained 
from nN processes. Their values were disputed by Entem and Machleidt [Til l30l|. and in this work we also try to 
investigate this question. 

An important point here is that we do not try to address details of the non-perturbative aspect of the NN 
interaction; instead, we deal only with the two-nucleon irreducible diagrams up to order g 4 , which should be identified 
with the irreducible kernel in a Bcthe-Salpeter equation, or the potential in a Lippmann-Schwinger equation. As 
suggested by Weinberg [3lJ, the usual power counting of ChPT should be applied to it. The reducible diagrams are 
afflicted with infrared divergences (or, in the language of old-fashioned time-ordered perturbation theory, small energy 
denominators) and, depending on the power counting adopted, a rcsumation of certain classes of diagrams have to 
be performed to all orders. The non-perturbative problem is dealt in the more generic framework of effective field 
theory (EFT for reviews see, for instance, Refs. j^M EE IBjl ) an d, despite of the progress it has achieved in the 
last years EHH |H E El E3 there are still controversies aboutpower counting schemes |3ll III III II^ . 

This work follows the same spirit as the ones from Kaiser et.al. [481 l49l |50| . Entem and Machleidt [l7]|, and Ballot, 
Rocha, and Robilotta [Ml , in the sense that we look only at peripheral waves (I > 3) where the Born approximation 
is expected to be valid |38t f5^) . With this procedure we try to investigate how far we are from the dream of establish, 
as in the case of the one-pion exchange (OPE) in the mid-60s, a mathematical structure for the two-pion exchange 
(TPE) which can be called "unique" , dictaded by chiral symmetry (by unique here we do not consider ambiguities 
arising from off-shell effects [k| or choices of field variables |37| , as discussed in Ref. [54] ) . 

We organize this paper as follows. In Sec. [H] we present the motivation of our previous works with some detail, 
showing the contrast between RB and HB calculations for the TPEP. There we also address the differences we found 
when comparing our expanding results with the HB expressions, and their origin are better (but still not completely) 
understood. An interesting outcome of this work is a quantitative measure of how different ranges of the potential 
contributes to the phase shifts, and this was made possible by the use of the phase function (or variable phase) 
method, described in Sec. IIIII This study, as well as the dependence of the phase shifts with some of the existing 
values for the LECs, is shown in Sec. II VI where we found some restrictions in order to achieve consistency with the 
experimental data. We close this work with our conclusions and remarks in Sec. 



II. RELATIVISTIC VERSUS HEAVY BARYON TPEP 

In our recent works 0, 0] we argued that the problem of convergence that appears in the HB formulation 
manifests itself in the TPEP at large distances. This was ilustratcd by comparing our basic loop functions with their 
1/m expansion, the latter not being able to reproduce the correct asymptotic behavior. In the end of this section we 
show that the same occurs to the TPEP decomposed in peripheral partial waves (J > 3). Before that, we present a 
discussion about the differences we found between our expanded results and the ones obtained in the HB framework 
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|171l48U49t 15(3 ]. It is worth to emphasize that they must in principle coincide, therefore, such discrepancies should not 
exist. To investigate them we revised some of our calculations, and updated results are presented in Sees. Ill 51 and 
III CI In order to fix the notation and help understanding the origin of these discrepancies we outline the main steps 
on the derivation of our TPEP. For further details, we direct the reader to our previous works Il5lll6l . For numerical 
evaluations, in this section we use the same values for the LECs used by Entem and Machleidt [TJ (last column of 
table EB. 
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FIG. 1: The irreducible two-pion exchange amplitude. 

The TPE contribution to the NN interaction can be parametrized in terms of ttN scattering sub-amplitudes, as 
ilustrated by the first graph on the right hand side of Fig. ^ O ne can describes the irN scattering in terms of two 
invariant amplitudes D ± (i/, t) and B^(v, t), 
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where m is the nuclcon mass, with initial and final momenta p and p' . respectively, while v = [(p + k) 2 — (p— k') 2 ]/4m 
and t=(k'— k) 2 are the usual Mandelstam variables. This allows the TPE amplitude to be written as 
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where the symbol /[■■■] represents the four-dimensional integration with two pion propagators, 

d 4 Q 1 1 

(2tt) 4 (fc 2 - [i 2 ) (k 12 - p 2 ) ' 



(2.3) 
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fi is the pion mass, Q = (k + k')/2, and the profile functions X's are loop integrals written in terms of the amplitudes 
D ± and B ± , 
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Using the variables 
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W = pi + p 2 = p'x 
we decompose the Lorentz structure of these profile functions as 



*= g [(Pi+K)-(P2+^)] 
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(2.6) 
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The evaluation of Ea. (|2.5|) to leading order gives rise to bubble, triangle, crossed box, and plannar box diagrams, the 
latter one containing an iteration of the one pion exchange (OPE). In order to obtain the irreducible TPE amplitude, 
this contribution has to be subtracted. Schematically it is represented by Fig. [3 where the symbol "+" on the second 
graph means that the two-nucleon propagator contains only its positive-energy projection. This projection is not 
uniquely defined 0,113 an d several prescriptions, like the Blankenbecler-Sugar y3(|, Equal-Time [57|. or the Gross 
|58| equations, could be used. Therefore, when comparing our expanded results with the HB ones, we should keep in 
mind that part of the differences might come from the particular prescription adopted. 
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FIG. 2: The plannar box diagram (left) and the positive-energy projection of the iterated OPE (right), which has to be 
subtracted to obtain the irreducible TPE amplitude. 

The iterated OPE can be recasted in the same structure as Ea. l|2.Hfl . allowing us to define the functions I's as 
the ones in Ea. (|2.5(l subtracted by the iterated OPE amplitude. The components of the potential, up to 0(q 4 ), can 
therefore be written in terms of these functions as 
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The dynamics of our rclativistic TPEP is classified in terms of three families of diagrams, according to their 
topology. The first line of diagrams in Fig. corresponds to the irreducible one loop graphs with vertices from the 
0(q 1 ) irN chiral Lagrangian, C^ N , with coupling constants at their physical values (family I). The second line (family 
II) contains two loop diagrams with an intermediate tttt scattering, while the third line (family III) comprises one loop 

(2) (3) 

graphs with vertices from C^-^ and C^n, &s well as other two loop graphs, parametrized in terms of irN subthreshold 
coefficients [Hflj . 





FIG. 3: Dynamics of the relativistic TPEP. The small black dots represent vertices from C^ N , the big shaded dots, the hit 
scattering amplitude, and the big black dots, the 7riV subthreshold coefficients. The latter contains other two loop contributions, 
as well as vertices from £^ and C^ N . 



B. revision of our two loop diagrams 



The expressions for our previous profile functions are given in appendices D-F of Ref. |l5j. While investigating 
the origin of the differences between our expanded and the HB results, we realize a small mistake in our two loop 
contributions (family II) . We checked that the revised calculation is independent on the parametrization of the pion 
field and leads to 
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where §a and f T are respectively the nucleon axial coupling and the pion decay constant, and IT, II t , II x , IT, IT, and 
II a are our basic loop integrals defined in Ref. These modifications alter the expressions of just two components 
of our non-expanded, relativistic potential: in Vq , the non-local (nonstatic) term 
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does not exist, while in V LS there are no contributions from family II at all. 
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C. 1/m expansion and analysis of discrepancies 



The second revision we made concerns the HB expansion of our loop functions. We found that the plannar and 
crossed box integrals (lit, and n x , respectively) which, in Ref. ^3], were given by 
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Therefore, the comparison of our expanded results with the HB expressions (Sec. X of Ref. mus t be updated 
by the following: 
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where L{q) is the usual logaritmic function, 
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and, from now on, q = \q\. 

We can notice that the number of different terms dropped from 14 in Ref. 0] down to 9. Now we have a better 
understanding about their origins: those marked with [• • -]t (six terms) are related to our choice of the Blankenbecler- 
Sugar [56j prescription to deal with the iterated OPEP, which is the same adopted in Refs. [II EE El il 111 . Diff erent 
prescriptions lead to different expressions only for these terms, as mentioned previously. In the HB calculations of 
Refs. l4Hl5(l . the plannar box diagram is expanded in powers of 1/m and then the iterated OPEP (given by Eq.(20) 
of Ref. [43) is identified and subtracted. A more detailed study of these aspects will be given elsewhere roll. 

The remaining three different terms, marked with [• • •]*, are originated from two loop contributions of |49j, part 
of them parametrized in our calculations through the subthreshold coefficients (family III). The one in Vr (and, 

consequently, Vg), Ea. H2.24[) . can roughly be identified with the unintegrated term Vj. of Ref. 01 > which was not 
considered while performing the above comparison. In order to show that this can be the case, notice first that their 
structures are similar. The HB expressions from Entem and Machleidt read |8a | 
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where n = \J a 2 /4 — /i 2 . Our equivalent terms are 
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In the first and second columns of table[I]we display, respectively, the values of Eas. (|2.35|l and l|2.33|l in configuration 
space, as a function of the internucleon distance r. We can see that the difference between them decreases with r, 
being around 20% at lfm, 10% at 2fm, and remaining in a few % for r > 3fm. As both expressions are numerically 
consistent at large distances one can say that this particular discrepancy is now understood. Looking at the third 
column, which shows the contribution of the term V^ a \ one can notice that such a difference is tiny compared to 
the whole two loop contributions, close to 3% at lfm and less than 1% for r > 2fm. To complete this analysis, we 
present in the fourth column the total value of the isoscalar tensor potential, which clearly shows that the whole two 
loop contributions (parametrized in our calculations inside family III) becomes less important at large distances. This 
behavior is already known, as ilustrated in Fig. 5(c) of Ref. 0|. Therefore the difference of using either Ea. (|2.33|l or 
Ea. (|2.35|l compared to the total Vr is numerically insignificant. 

The two remaining discrepancies, which appear in Wc, are more difficult to trace. Analytically the identification 
of terms is not as clear as in Vt and, in principle, does not appear to be equivalent. Despite of that, we follow the 
same numerical procedure as before, identifying W^? in our calculations as 



W, 



(b) 



L(q) 



(15 + 7ffi) 



18432tt 4 /^ 



(V + <7 2 ) 



2p 2 (2 M 5 



+ g 2 )-3/5( 5 2 -l)(4 M 2 + q 2 ) 
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TABLE I: Values of the isoscalar tensor potential (MeV) in configuration space. The first, second, and third columns are 
defined and explained in the text, while the last column is its total value. 



r( fml 


V^ b \ Ea.j2.SSI 


V^ b) , Ea.j2.SSt 


vi a) 


Vt (total) 


1 .00 


-4.5479 E+00 


-5.4716 E+00 


-2.4737 E+01 


-3.9311 E+01 


2.00 


-2.2078 E-02 


-2.4062 E-02 


-1.2008 E-01 


-3.2551 E-01 


3.00 


-6.6798 E-04 


-7.0387 E-04 


-3.6332 E-03 


-1.5170 E-02 


4.00 


-4.0679 E-05 


-4.2189 E-05 


-2.2126 E-04 


-1.2858 E-03 


5.00 


-3.5582 E-06 


-3.6576 E-06 


-1.9353 E-05 


-1.4573 E-04 


6.00 


-3.8725 E-07 


-3.9586 E-07 


-2.1063 E-06 


-1.9562 E-05 


7.00 


-4.8725 E-08 


-4.9619 E-08 


-2.6502 E-07 


-2.9308 E-06 


8.00 


-6.7941 E-09 


-6.8999 E-09 


-3.6954 E-08 


-4.7418 E-07 


9.00 


-1.0225 E-09 


-1.0363 E-09 


-5.5614 E-09 


-8.1184 E-08 


fO.OO 


-1.6319 E-10 


-1.6514 E-10 


-8.8763 E-10 


-1.4515 E-08 



d>g\{2n 2 +q 2 ) 



(9 2 A 



1) (4/z 2 



while the HB expressions are given by 



w, 



2q 6 f°° , ImWk b) (ia) 

da ^ — - 



3a( 



2jU a 5 (a 2 + q 2 ) ' 

^-3 f dx [g\{2f a 2 ) + 2{g\ 1)kV] 
^Jtv) Jo 



x | — 3n 2 x 2 + Qkx^J /i 2 + n 2 x 2 In 
+g\(a 2 - 2k 2 x 2 - 2^i 2 ) 



KX + yj /.I 2 + K 2 X 2 



/' 



6 K 2 X 2 \ K 2 J 



3/2 



(2.36) 



(2.37) 



In 



KX + \J (I 2 + K 2 X 2 



(2.38) 



In table UTI we show the analogous of table Q] for Wc, where one cannot observe the same consistency as before: at 
large distances (r > 4fm), our result for (what we initially supposed to be) is almost twice larger than the HB 

one. Fortunately this difference is not numerically important for the total value of Wc (as one can check by looking 

at the last column), but is somehow relevant for the whole two loop contributions (Wq + W@). This numerical 
comparison indicates that these two mentioned discrepancies still persist, and should be considered the last small 
detail to be solved in order to reach an unique structure for the 0(q 4 ), chiral two-pion exchange component of the 
nucleon-nucleon potential. 



D. clean comparison between relativistic and HB TPEP 



Here we present the consequences of expanding our TPEP in powers of 1/m. In order to isolate only this particular 
effect, we use the fact that our expanded results have to coincide with the HB expressions. As the discrepancies in 
Eas. H2.23|) - (|2.3U|) are not actually related to the problem of the HB expansion, and using the fact that they can easily 
be identified in our non-expanded expressions, we simply ignore them in our relativistic calculations from now on. 
With this procedure we can also compare our phase shifts computed in configuration space with the ones from Entcm 
and Machleidt performed in momentum space. 

In Fig. 0] we plot the ratios of the HB over the RB components of the TPEP. Even though the distance considered 
is quite large for hadronic interactions, we will show latter in Sec. IIIII that peripheral phase shifts at sufficiently small 
energies are sensitive to distances as large as lOfm. One can see that the pattern observed in the 1/m expansion of 
our loop functions 0] reflects more on the isovector (-) channels, but produces deviations smaller than 10% on the 
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TABLE II: Values of the isovector central potential (MeV) in configuration space, analogously to what was done in table[l]for 
the isoscalar tensor component. 



r( fml 






yv c 


Wc (total) 


1.00 


-5.8826 E+00 


-5.5372 E+00 


-7.6350 E+01 


-5.3305 E+01 


2.00 


-3.3949 E-02 


-2.4611 E-02 


-3.2322 E-01 


-4.5217 E-01 


3.00 


-1.2397 E-03 


-8.0784 E-04 


-8.9086 E-03 


-2.7783 E-02 


4.00 


-9.0177 E-05 


-5.5440 E-05 


-5.0169 E-04 


-3.0848 E-03 


5.00 


-9.2668 E-06 


-5.4899 E-06 


-4.0838 E-05 


-4.4010 E-04 


6.00 


-1.1654 E-06 


-6.7289 E-07 


-4.1468 E-06 


-7.1764 E-05 


7.00 


-1.6697 E-07 


-9.4596 E-08 


-4.8701 E-07 


-1.2702 E-05 


8.00 


-2.6176 E-08 


-1.4617 E-08 


-6.3325 E-08 


-2.3757 E-06 


9.00 


-4.3816 E-09 


-2.4189 E-09 


-8.8692 E-09 


-4.6224 E-07 


10.00 


-7.7072 E-10 


-4.2156 E-10 


-1.3135 E-09 


-9.2646 E-08 



1.5 




r(fm) 



FIG. 4: Ratios between the HB and RB components of the TPEP, in configuration space. 



isoscalar (+) ones. Note that in the isoscalar central potential it remains below 5% up to lOfm. Wc omit the ratios 
of the spin-spin components, which arc almost the same as for the corresponding tensor ones. Before commenting 
further on these results, it is interesting to look how they manifest in the partial wave decomposition. 

Figs. [SJ El andJTJshow, for F, G, and H partial wave channels, respectively, the values of the RB and HB versions of 
the TPEP multiplied by a common factor exp(2^r), as well as their corresponding ratios. The notation is similar to 
the one used for phase shifts — for a given total and orbital angular momenta J and L, the single channel potentials 
are identified by J = L waves, while for coupled channels the 2x2 matrix potential has diagonal components 
[ 3 (J — 1),/, 3 ( J + l)j] and non-diagonal component ej. A chiral TPEP in configuration space has the drawback of 
being divergent at the origin and, in order to avoid numerical complications while computing phase shifts, wc multiply 
it by a phenomenological regulator, 



V 



TPEP, reg 



1 



4 



l^PEP, div j (2.39) 



which is similar to the one used by the Urbana |62| and Argonne |63j, |6J] potentials. For the regulator we adopted 
the value from AV14 c = 2.0fm -2 . We will address this question with more details in Sec. Ivl 

Besides the non-diagonal potentials, the ratios clearly show the failure of the HB in describing the correct asymptotic 
behavior of the TPEP. Such deviations get more dramatic for 1 F 3 , 3 Gs, 3 G*4, 3 G5, and 1 H 5 waves (not by accident), 
being significant already at r = 5fm. From Fig. 0] it seems strange that differences of the order of 35% turn into 
~ 50 — 100% observed for these waves. The reason is that they are channels with total isospin T — 0, leading to 
cancelations between isoscalar and isovector components. As an example, for total spin S = states the tensor and 
spin-orbit components don't contribute and the potential is given by 
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V [ l Lj] = V+ + (4T - 3) V c - 3 [V+ + (4T - 3) V^] = V + - 3 V c - 3 [V+ - 3 V^] . 
At r = 6fm one has, for the RB and HB cases, 



(2.40) 



V \^Lj\ RB = [ - 150.0 - 3(-78.7)] - 3 [30.1 - 3(4.28)] w 86.1 - 51.8 w 34.3 
v [ 1l j]hb = [- 150.7 -3(-88.8)] - 3 [30.4 - 3(3.78)] w 115.7 - 57.2 w 58.5 



(2.41) 
(2.42) 



(in units of 10 _6 MeV). From Eqs.lEHU) and lt2~4*2l one can understand the mechanism that leads to a deviation 
of roughly 70% in the final result: it comes first from the structure [(isoscalar)-3(isovector)] (which generates more 
than 30% of deviation from the central component), and is further increased by the sum of central and spin-spin 
contributions. 
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FIG. 5: Values of the RB and HB potentials [multiplied by exp(2fj,r)], as well as their ratios, for F waves. 
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FIG. 6: Same as in Fig. for G 
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FIG. 7: Same as in Fig.0 for H waves. 
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III. THE PHASE FUNCTION METHOD 



In the calculation of phase shifts we employ the formalism of the variable phase, well studied by a couple of authors 
in the early sixties EE EE E3, EE EE EE but, surprisingly, not very popular in the nuclear ph ysic s community (in the 
last decade, from the best of our knowledge, only a few works took advantage of this method [42l IHH IH2I Flj ). Due to 
its simplicity, here we outline the main steps to obtain the phase function equation, where its physical interpretation 
immediately arises. A more detailed description, including scattering with tensor forces, as well as extensions of these 
ideas to other processes, observables, and kinematics can be found in a review paper by Babikov |68| . as well as in 
the book from Calogero EE ■ 

Here we consider the case of a central potential with a regular behavior, i. e., the ones which behave as r~ n close 
to the origin, with n < 2. The radial Scrodinger equation is written as 

<(r)+ [k 2 -L(L + l)/r 2 -U(r)] u L {r) =0, (3.1) 

with U(r) = 2m rec iV(r), and k = yJ2m re d£, V(r) being the central potential, £, the energy eigenvalue in the 
Schodinger equation, m rec i, the reduced mass of the system, and ul{t), the radial wave function. 

The basic idea of the phase function method consists in introducing the functions 5l{t) and Al(?~) in such a way 
that the radial wave function and its derivative are written as 

u L {r) = A L (r)[cos5 L (r)j L (kr) - sin 5 L (r) n L (kr)] , (3.2) 
u' L {r) = A L {r)[cos8 L {r)$ L {kr)-sm8 L {r)h' L {kr)} , (3.3) 

where Jl{x) and fiL(x) are the usual spherical Bessel functions, Jl(x) and iil(x), multiplied by the argument (so called 
Riccati-Besscl functions). Ea. (l3.2|) is just a general paramctrization of wl(^), while Ea. l|3.3|l imposes the condition 

A' L (r) [cos 5 j L(kr) — sin SL{r) hL(kr)] — S' L (r) Al(t) [sin<5i(r) jh{kr) + cos5l(t) fiL(kr)] =0. (3.4) 

Substituting Ea. (|3.2() in Ea. l|3.1[) . using the condition l|3.4ll . and trigonometric relations we end up with a first order 
diffential equation for 8l{t) 1 the "phase function": 



S' L (r) 



U(r) 



cos 



8l (r) j L (kr) - sin S L (r) n L (kr) 



(3.5) 



Notice that beyond the point where it is assumed to be the range of the potential this function becomes a constant, 
and the form of Ea. 13.2fl tells us that its value is precisely the phase shift. Furthermore, the condition l|3.4|l assures 
the continuity of the derivative of the wave function at this point. But the phase function means more than that — 
from Ea. (|3.5|l one can easily check that, for a given distance R, the quantity 8l(R) indeed is the phase shift associated 
to a potential cut at this point, V(r) 0(R — r). 




FIG. 8: Projections of the AV14 NN potential into 
So, D2, and 1 G f 4 partial waves. 



FIG. 9: Normalized phase function at Ei a b — 50MeV 
associated to the partial waves from Fig. |H] 



13 



This property is particularly interesting to assess the importance of different regions of the potential. For example, 
in Fig. E|we show the normalized phase function, 6l(t)/6l(oo), for nucleon-nucleon scattering in 1 Sq, 1 I?2, and 1 C?4 
partial waves, with the laboratory energy at 50MeV. The corresponding potentials, displayed in Fig. 03 are partial 
wave projections of the Argonne V14 potential |6^|. Note that, for 1 Sq, the repulsive potential in the region r < 0.8fm 
is described by a decrease in the value of the phase function up to a minimum at this point. Beyond that the potential 
is attractive, resulting in positive increments to 5l{t). For the other waves, the phase function also shows nicely the 
effect of the centrifugal barrier which, for sufficiently large L and small energy, supress the dependence of the phase 
on the short range structure of the potential. In the next section we take advantage of this feature to determine the 
ranges of relevance when extracting the phase shifts. 

IV. LECS AND PHASE SHIFTS 

In this section we present our studies for peripheral phase shifts. The calculations are performed in configuration 
space using the phase function method described in the previous section, which is fully equivalent to the Schrodinger 
equation. The final NN potential is the sum of TPEP and OPEP and we disregard the exchange of three or more 
pions, which are much more weaker and smaller in range [3^. F5 , l73[ . We use the Fourier transform of the same 
charge-dependent OPE for neutron-proton scattering used by Entem and Machleidt [l7j . 



Vope = -KGM + (-i) t+1 2K(m), 

with n„o = 134.9766 MeV, fx^± = 139.5702 MeV, and the function V n (ji) defined as 



(4.1) 



4/2 



i + A 



(H 5 



Sia + <r (1) • <x (2) 



(4.2) 



where Si 2 is the usual tensor operator in configuration space. This is the same charge-dependent OPEP used in 
Refs. |6J,|7J|, provided we set the scaling mass m s and the couplings as 



fpp — fnn 



1 



(4.3) 



For our TPEP we use the average nucleon and pion mass from the AV18 potential [64| . m = 938.9190 MeV and 
H = 138.0363 MeV. The remaining constants besides the LECs are he = 197.3270 MeV.fm, f n = 92.4 MeV, and 
qa = 1.3187, the latter being higher than the experimental value in order to accommodate the Goldberger-Treiman 
discrepancy. In this work we did not consider the quadratic spin-orbit term, given by Eq.(22SJ. Its effect on the 
phase shifts is expected to be insignificant for the waves we are considering as we could estimate, for instance, by 
switching on and off the L 2 and (L ■ S) 2 components of the Argonne V14 potential and comparing the results. 



TABLE III: Values for the dimension two (d) and three (di) LECs considered in this work, given respectively in GeV 1 and 
GeV" 2 . 



LEC 


Ref. [21j 


Ref. [23J (fit 1) 


Ref. [24] 


Ref. [22] 


Ref. [1JJ 


Cl 


-0.94 


-1.23 


-0.81 


-0.76 


-0.81 


C2 


3.20 


3.28 


8.43 


3.20 


3.28 


C3 


-5.40 


-5.94 


-4.70 


-4.78 


-3.40 


C4 


3.47 


3.47 


3.40 


3.96 


3.40 


di + d 2 


2.40 


3.06 


3.06 


2.40 


3.06 


d 3 


-2.80 


-3.27 


-3.27 


-2.80 


-3.27 


d 5 


1.40 


0.45 


0.45 


1.40 


0.45 


di4 — di5 


-6.10 


-5.65 


-5.65 


-6.10 


-5.65 



Concerning the LECs, in table lllll we display their values from five different studies. The first column was extracted 
from a work from Mojzis plj where he presented the first complete 0(q 3 ) expressions for ttN scattering amplitude and 
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pinned down the LECs using the Karlsruhe-Helsinki [26J threshold parameters in combination with the pion-nucleon 
a term. The second column was determined by Fettes et.al. through a direct fit to S- and P- wave ttN phase 
shifts between 40 and 100 MeV pion lab momentum, usin g th e Karlsruhe partial wave analysis [22j (fit 1). The 
third column shows the results from Buttiker and Meifiner j24|, where the irN amplitude was reconstructed inside 
the Mandelstam triangle via dispersion relations, then the LECs were obtained through a fit to this amplitude, in 
two different points (because their results for the dimension three LECs were not consistent, we adopt here the same 
values from the second column). These three works rely on old Karlsruhe dispersion analyses, which have presumably 
some inconsistent data included l75j |87l. Different from them, the dimension two LECs on the fourth column were 
determined by the Nijmegen group [29j from their new partial wave analyses of 5109 proton-proton and 4786 neutron- 
proton scattering data below 500 MeV. As the dimension three LECs were not considered in their work, for this set we 
adopted the same values from Mojzis. Finally, on the last column we have the values used by Entem and Machleidt 
in their calculations of neutron-proton phase shifts [T^ . 

In the calculation of phase shifts we use the familiar nuclear-bar convention |78[ , where the coupled S matrix with 
total angular momentum J is parametrized as 



a _{ e 2 "- 7 " 1 cos2ej i e MSj-i+s J+ i) sin 2ej \ 

bj ~ y ie HSj-i+6j+i) S in2ej e 2 ^ 1 cos2e, 7 J' (4 ' 4) 

Before addressing their dependence on different sets of LECs it is interesting to take advantage of the phase function 
formalism fSec, IIIl|l to identify the most important ranges of the potential in the determination of these quantities. In 
Figs. El ^2 and 1 121 we plot the distances where the phase function 6l{t) reaches 10%, 20%, and so on, up to 90%, of 
the final value, <5l(oo), as functions of the laboratory energy, -Elab- As in Sec.[n] we use here a potential with LECs 
from Entem and Machleidt (last column of table IIIl|l and similar curves are obtained for the other sets. These figures 
show, quantitatively, the sensitivity of phase shifts on large distances at sufficiently low energies. We also notice that 
the range, delimited by the points where the phase function reaches 10% and 90% of the final value, gets narrower 
and closer to the inner region of the potential as we go higher in energies, but not smaller than lfm (1.5fm) for F (G 
and H) waves, a clear manifestation of the centrifugal barrier. 





FIG. 10: Distances, for F waves, where the phase shift (or mixing parameter) achieves X% of its value, as functions of the 
laboratory energy. The curves correspond, from bottom to top, to X —10, 20, • ■ ■, 90, respectively. 
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FIG. 11: Same as Fig. GUI for G waves. 
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FIG. 12: Same as Fig. [T@ for H waves. 

Figs. El EH an d show, respectively, the phase shifts associated to F, G, and H waves, as functions of the 
laboratory energy. The dotted lines correspond to LECs from Mojzis (column 1), dot-dashed, from Fettes et.al. 
(column 2), dashed, from Biittiker and Meifiner (column 3), solid, from the Nijmegen group (column 4), and thick 
solid, from Entem and Machleidt (column 5). The dark and light lines represent, respectively, calculations using the 
RB and HB formalism. As it is clear, from NN scattering phase shifts one cannot observe any significant difference 
between RB and HB results — despite the contrast between them, shown in Sec. Ill Dl in NN scattering one also has 
to add the OPEP, which gives a large contribution and makes such discrepancies harder to observe. Exeptions are 
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the F wave results using the LECs from Fcttes et.al, but their predictions are not very close from what one should 
expect. 

It is noticeable the large variations of the phase shifts with different sets of LECs, predominantl y d ue to C3 (the 
leading order term of the nucleon axial polarizability |x3 ) ■ As pointed out by Entem and Machleidt |l7j , consistency 
with the experimental analysis fFigs. fTTA 1171 and llSj l favors a smaller value for this particular LEC, C3 ~ —3.4. We 
will discuss this issue later in Sec. 




50 100 150 200 250 300 50 100 150 200 250 300 

E^MeV) E^MeV) 
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FIG. 13: Phase shifts and mixing parameter for F waves, as functions of the laboratory energy. The dotted, dot-dashed, 
dashed, solid, and thick-solid lines correspond to LECs from table ITTTI columns 1, 2, 3, 4, and 5, respectively. The dark curves 
are results from the RB formalism, while the light ones, from the HB expressions. 
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FIG. 14: Phase shifts and mixing parameter for G waves, as functions of the laboratory energy. Notation is the same used in 
Fig. ESI 
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FIG. 15: Phase shifts and mixing parameter for H waves, as functions of the laboratory energy. Notation is the same used in 
Fig. ESI 
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In Figs. ^3 El and 1181 wc compare the phase shifts associated to the LECs from column 4 (dark, dashed curves) 
and column 5 (dark, solid lines) of table ITTTI with the Nijmegen partial wave analysis (circled line) |79j |. It is important 
to emphasize that, in Refs. jlSESli the Nijmegen group employ the TPEP up to 0(q 3 ) to determine the dimension 
two LECs ci, C3, and C4. For this reason, we also plot the phase shifts using the TPEP up to this order with their 
LECs, indicated by the (light) dot-dashed curves. 
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FIG. 16: Phase shifts and mixing parameter for F waves, as functions of the laboratory energy, compared with the Nijmegen 
partial wave analysis (circled line) |79|. The solid and dashed (dark) lines correspond to OPE plus TPE up to 0(q 4 ) using, 
respectively, the LECs from Nijmegen (column 4) and Entem and Machleidt (column 5), while the dot-dashed (light) curves 
are results from OPE+TPE up to 0(g 3 ), using the LECs from Nijmegen. 
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FIG. 17: Phase shifts and mixing parameter for G waves, as a function of the laboratory energy. Notation is the same used in 
Fig.EO 



21 



50 100 150 200 250 300 50 100 150 200 250 300 

E^MeV) E^MeV) 




50 100 150 200 250 300 50 100 150 200 250 300 

E^MeV) E^MeV) 




50 100 150 200 250 300 



E UB (MeV) 



FIG. 18: Phase shifts and mixing parameter for H waves (as a function of the laboratory energy). Notation is the same used 
in Fig. OH 
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V. COMMENTS AND SUMMARY 



As mentioned in Sec. Ill Dl the calculation of phase shifts requires the regularization of the potential near the origin 
and, as a consequence, makes our calculation dependent on the value of the regulator. This procedure is analogous to, 
in momentum space, introduce a cutoff in the dynamical equation, and seems to be unavoidable when dealing with two 
or more nucleoli systems. There are, though, attempts to reformulate the problem using dimensional regularization 
j8f| or variants |l2l H3 | in the literature. One should emphasize that even if the chiral potential is finite in momentum 
space, its representation in configuration space shows a divergence at the origin of the type r~", where n is some 
integer. It simply maps the fact that a chiral potential has a limited region of applicability, increasing as a polinomial 
in q (in contrast, for instance, with phenomenological one-boson exchange potentials). The advantage to work in 
configuration space is that statements about ranges and nucleon sizes can be made more transparently, and also it is 
more suitable to be extended to few bodycalculations for nuclei up to A = 10, e.g., by means of the Variational or 
Green Function Monte Carlo techniques [81j . 

We introduced a regulator given by Ea. (|2.39|l . with c = 2.0fm~ 2 , in our calculations. In principle there is not a 
real compelling justification to use such a value. With higher values, however, the phase shifts using the LECs from 
Fettes et.al., or from Mojzis, become even more unrealistic. For lower values the situation interchanges, as ilustrated 
in Fig. ^1 but at the expense of large modifications in the strength of the potential. For example, using c = l.Ofm -2 
one has a reduction of almost half of the "bare" value of the potential at r = 1.4fm. With c = 2.0fm -2 at the same 
point, it reduces roughly 8%. From the same figure we see that, even with this regulator, the sets with larger absolute 
values of C3 are still far from PWA results fFigs. IT?)lll8(l . There might be two explanations for this to happen. One, 
that F waves are not peripheric enough to avoid the dependence on short distances of the potential, but we check 
that similar behavior holds for 1 G*4 wave as well. Also, one can rule out this option by calculating the phase shifts in 
Born approximation, as shown in Fig. 121)1 where there is no need to regulate the potential near the origin. 





100 150 200 250 300 

E UB (MeV) 



FIG. 19: 3 F 4 phase shifts, using c = l.Ofm" 2 for the 
regulator. Notation is the same used in Fig. 1131 



FIG. 20: Born approximation for F4 phase shifts. No- 
tation is the same used in Fig. 1131 



The other possible explanation rely on EFT ideas to deal with divergences in the dynamical equation, implemented by 
the introduction of a cutoff. To be consistent, the same cutoff should also regularize loop integrals for the irreducible 
two- nucleon diagrams that constitute the kernel (or the potential) of the dynamical equation [88j . This procedure can 
weaken the contribution of pion loops at short and intermediate distances, enough to allow larger absolute values for 
C3 |38|. We do not address this alternative in the present work, as we keep the possibility of renormalize the dynamical 
equation by other means, as the ones mentioned before. 

A further remark on the regulator is that our calculations, using the same LECs from Entem and Machleidt, are not 
very sensitive to its variations in the range l.Ofm -2 < c < 3.0fm -2 , except for 3 F"3 and 3 F4 waves. Using c = 2.5fm -2 , 
we have all of our results very close to theirs 0| . 

One comment on Figs. I16H18I is in order. We can see that using the LECs from the Nijmegen group on the 
expressions of the 0(g 4 ) TPEP does not give a good description of results from PWA, but this procedure does not 
make sense, as these LECs were determined using a potential to 0(q 3 ) |2^,|2^|. It is well known that the LECs pinned 
down by any particular process depend on the order one is working on. Formally the difference between using an 
0(q 3 ) and 0(q 4 ) TPEP should be of higher order, but in practice it is sizeable (see, for instance, the ttN scattering 
analysis from Ref. J8j| ). We can check that the too strong attraction, which arises if one uses the LECs from the 
Nijmegen group in the 0(q A ) potential (Nij 4 ), is actually not that attractive if one considers it up to 0(q 3 ) (Nij 3 ), 
bringing their phase shifts closer to the ones from PWA. Fig. 1211 ilustrates this point for 1 i 7 3 wave. The thick-solid 
curve is the potential using the LECs from Entem and Machleidt, while the dashed and solid curves are the Nij 3 
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and Nij 4 potentials, respectively. For comparison purposes, we also include the OPEP (light, dotted curve) and the 
phenomenological Reid 93 potential (dot-dashed curve). 
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FIG. 21: F3 potential. Notation is given in the text. 



From this figure and what we discussed above one can say that the argument given by Entem and Machleidt |l7ll30l| . 
that the LEC C3 from the Nijmegen group leads to a much stronger potential as a result of a cutoff in the potential at 
1.4fm or 1.6fm (and set to zero until the origin) could have some meaning only i/such a value were extracted from an 
0{q A ) TPEP. Evidently this is not the case, primarily because the energy-dependent boundary conditions of Nijmegen 
PWA at r = 1.6fm does not mean that the potential is zero inside this range. Besides that, the Nij 3 potential shows a 
repulsion beyond r ~ 1.2fm, as can be seen in Fig.[^] The (now, not too large) disagreement of Nij 3 must have other 
explanations. From our point of view, the reason is that the LECs are parameters in the TPEP beyond 1.6fm which 
are determined when the fit to scattering data gives the best x 2 with a minimum number of boundary conditions. 
This involves sums also over lower partial waves, where the nonperturbative aspect of the nucleon-nucleon interaction 
is more pronounced. Therefore, the way short distances are parametrized can contaminate the actual value of these 
LECs, contributing as a systematic error. Moreover, contact terms can start contributing at 0(q 5 ) or 0(q 6 ) to F 
waves, which can change the the behavior of phase shifts even at not so large energies j85j. When fitting to data, 
these terms are inherently taken into account. 

To summarize, in this work we address the discrepancies we found when comparing our 1/m expanded results with 
the expressions from the HB formalism. It forced us to revise not only our 1/m expansions, but also our two loop 
diagrams, and now we have a better understanding about the origin of the remaining discrepant terms. There arc, 
however, two terms in the central isovector component, originated from two loop calculations, which could not be 
understood. It will be necessary to solve them if claims for higher order calculations of the chiral TVTV potential 
becomes imminent. In a second stage we calculate the phase shifts for peripheral waves (F, G, and H), where short 
distances are supposed to play a minor role. This was shown in Figs. HUH 121 where the ranges relevant to phase 
shifts where determined using the method of the phase function. Our results for phase shifts reveals that, due to 
the overwhelming contribution of the OPEP, any difference between the RB- and HB- TPEP is barely noticed. It 
is possible that these differences might be relevant in processes where the OPEP is cancelled, and the dominant 
tail of the interaction, given by the TPEP, for instance, in low-energy peripheral scattering of nucleons off isoscalar 
targets [84|. A more relevant question, there are large variations of the phase shifts with the existing values of LECs, 
predominantly due to C3, which dictates the strength of the attraction of the potential. Consistency with PWA 
favors smaller values, in magnitude, than the existent ones in the literature. There are, however, conceptual details 
mentioned before which prevent, at the moment, to constrain these LECs with a satisfactory precision. They should 
be well understood and quantified if one aims to extract values with good precision from TV TV scattering. 
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